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This chapter describes the progress made during the past three decades in the 
finite size scaling analysis of the critical phenomena of the Anderson transition. 
The scaling theory of localisation and the Anderson model of localisation are 
briefly sketched. The finite size scaling method is described. Recent results for 
the critical exponents of the different symmetry classes are summarised. The im- 
portance of corrections to scaling are emphasised. A comparison with experiment 
is made, and a direction for future work is suggested. 

1.1. Introduction 

Originally, the phenomenon of localisation is a property of quantum mechanical 
wave functions bound in potential wells of finite range. At infinity, where the 
potential vanishes, the wave functions decay exponentially for negative energies 
indicating that the probability of finding the particle far from the potential well 
vanishes. This is called "potential localisation". It had already been suggested 
in the 1950s that potentials with infinite range can also support the existence of 
localised wave functions at positive energies provided that the spatial variation of 
the potential is random. This localisation phenomenon is due to destructive inter- 
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ference of randomly scattered partial waves and is now referred to as "Anderson 
localisation" . The most important physical consequence of Anderson localisation 
is the suppression of diffusion at zero temperature, which was conjectured by P.W. 
Anderson in his seminal paper. 1 Perhaps, the most striking example of Ander- 
son localisation is in one dimensional random potentials where all the states are 
localised, irrespective of their energy. The study of one dimensional localisation 
was pioneered by Mott and Twose 2 and by Gertsenshtein and Vasilev. 3 It can be 
treated exactly and has been the subject of several reviews. 

In higher dimensions the problem is more subtle, with the possibility of energy 
regions corresponding to localised states only, and to extended states only, sepa- 
rated by critical energies, called "mobility edges" . The zero temperature and zero 
frequency electrical conductivity <t of the system vanishes if the Fermi energy is 
located in a region of localised states. In the region of extended states, d ^ 0. In 
the absence of interactions, the system is an electrical insulator in the former case 
while in the latter case metallic conductivity is expected. It was conjectured in a 
seminal work 7 that this metal-insulator transition exists only in three dimensions, 
while in dimensions d < 2 systems are always insulating. This conjecture was based 
on the hypothesis of one parameter scaling of the conductance g(L) of a system of 
size L, i.e. that the dependence of the conductance on system size can be described 
by a beta-function, 

that depends only on the conductance. The behaviour of /3(g) with g was con- 
jectured based on perturbation theory in the limits of weak and strong disorder 
(large and small conductance), and assuming continuity and monotonicity in be- 
tween. Moreover, according to the scaling theory, at the mobility edge a continuous 
quantum phase transition between an insulator and a metal occurs accompanied 
by the power law behaviour of physical quantities, described by critical exponents, 
that is typical of critical phenomena at continuous phase transitions. The critical 
exponents of the conductivity 



<jo~(E-E ) s , (1.2) 



and the localisation length 



cT~(£ (1-3) 

were predicted to obey Wegner's previously conjectured scaling law 8 

s={d-2)v. (1.4) 

While this work was a great leap forward in our understanding of Anderson locali- 
sation, it remained to establish the validity of the central assumption of the theory, 
namely the one parameter scaling hypothesis. 
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This question was addressed numerically by simulating the Anderson model 1 of 
disordered quantum systems which consists of a delocalising kinetic energy mod- 
elled by a hopping term V and a localising random potential energy €j, commonly 
assuming a white noise distribution of width W , on a discrete square lattice {j}, 



where S denotes the nearest neighbours of the lattice site j. Such simulations allowed 
the one parameter scaling hypothesis to be verified with a reasonable numerical 
precision, 9-12 in the center of the band, at energy E = 0, and also to confirm the 
prediction s = v. 9 The critical disorder in three dimensions was initially found to be 
W C (E = 0) = 16 ± 0.5 while s = v = 1.2 ± 0.3. Although this latter value seemed to 
be consistent with v = 1 there were subsequently substantial doubts about whether 
or not this was indeed the case. It was found necessary to improve the precision of 
the estimate of the critical exponent and to study in detail and with high precision 
the conditions for the validity of the one parameter scaling hypothesis. Later, it was 
found that the exponent was in fact not unity and this intriguing discrepancy was 
the reason for numerous further numerical as well as analytical efforts, especially 
since the experimental situation was also far from clear. 13 

In the following sections, we briefly review the development of the finite size 
scaling analysis of the Anderson transition, paying particular attention to the role 
of symmetry and the estimation of the critical exponents. We stress the importance 
of the taking proper account of corrections to scaling, which has been found to be 
essential in order to estimate the critical exponents precisely. Finally, we tabulate 
the "state of the art" estimates for the critical exponents of the different universality 
classes. 

Some of the early results have been described in previous review articles. 14-17 
1.2. The Anderson Model of disordered systems 

In this section we briefly explain the Anderson model of localisation. We generalise 
Eq. (I1.5[) in order to describe more general systems with different symmetries. The 
most general form of Eq. (|1.5[) is 



The states | that are associated with the sites of a regular lattice j — usually 
for simplicity a square lattice is assumed — are assumed to form a complete set 
such that (jfi | j'fj,') = Sjj'Sfj,^' . Indices fi denote additional degrees of freedom 
associated with the lattice sites which lead to several states per site. If there are 
n states the above Hamiltonian describes Wegner's n-orbital model. 8,18 In general, 
the potential energies ej^ and the hopping integrals Vj^ji ^ are random variables 
described by some statistical distributions. 




(1.5) 




(1.6) 
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If the energy bands emerging due to the broadening by the kinetic terms are 
not strongly overlapping, we may use the single band approximation Eq. (|1.5|) . In 
addition, if we assume that sufficiently close to the Anderson critical point the 
critical phenomena are universal, i.e. independent of the microscopic details of the 
system, Eq. (|1.5|) is the simplest model that can describe the critical behaviour 
at the Anderson transition. In principle, these assumptions have to be verified a 
posteriori, and to some extent this has indeed been done during the past decades. 

If the Anderson transition is a genuine phase transition, the critical behaviour 
can be expected to depend only on symmetry and dimensionality. For a disordered 
system, spatial symmetry is absent and only two important symmetries remain: in- 
variance with respect to time reversal, and invariance with respect to spin rotations. 
Three symmetry classes are distinguished^: the orthogonal class which is invariant 
with respect to both time reversal and spin rotations, the symplectic class which is 
invariant with respect to time reversal but where spin rotation symmetry is broken, 
and the unitary class where time reversal symmetry is broken. Note that, if time 
reversal symmetry is broken, the system is classified as unitary irrespective of its 
invariance, or otherwise, under spin-rotations. 

When the kinetic energy parameter V is a real number, the Anderson model 
Eq. (| 1 . 5[) is time reversal invariant and belongs to the orthogonal class. In this 
case, universality has been verified by showing that a Gaussian, Cauchy and a box 
distribution of the disorder potential give the same critical exponents. 20 

When the kinetic terms Vjj/ become complex, the system is no longer time 
reversal invariant and thus belongs to the unitary class. This can be physically 
realised by applying a magnetic field. Then, the hopping term has to be replaced 
by the Peierls substitution 



where the vector potential A describes the magnetic field, B = VxA. Two different 
unitary models can be constructed using the Peierls Hamiltonian, namely a random 
phase model which is characterised by 



with the uncorrelated phases tpjji as random variables, and a model of a uniform 
magnetic field that leads to a similar expression for the kinetic term but with cor- 
related phases. Whether or not these two unitary models have the same critical 
behaviour has been the subject of numerous studies. 

In the presence of spin-orbit interaction, spin rotation symmetry is broken. The 
simplest Hamiltionian for such a symplectic case is 21-23 



a In fact, the classification is more complicated. 16 ' 17,19 However, for the present purposes, the 
following classification is sufficient. 




(1.7) 



Vjf = Vexp(iip jr ) , 



(1.8) 




(1.9) 
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where Uj ji is an SU(2) matrix. This model describes a two dimensional electron 
system in the presence of Rashba 24 and Dresselhaus 25 spin-orbit couplings. 

If the Anderson transition is a genuine quantum phase transition, we expect 
that the critical behaviour is universal and that the critical exponents depend only 
on the symmetry class and the dimensionality. 

1.3. Finite size scaling analysis of the Anderson transition 

In principle, phase transitions occur only in the thermodynamic limit, i.e. in an 
infinite system. In practice, computer simulations are limited to systems of small 
size. This necessitates an extrapolation to the thermodynamic limit. This extrapo- 
lation is far from trivial. It requires a numerically stable procedure which, at least 
in principle, allows control of the errors involved. This is especially the case when 
the goal is precise estimates of the critical exponents. Finite size scaling is such a 
procedure. 

1.3.1. Finite size scaling 

The raw data for the finite size scaling procedure is some appropriate physical 
quantity in a system of finite size. For some physical quantities it may be necessary 
to take a statistical average. An example is the two terminal conductance where 
an average over a large number of realisations of the random potential is required. 
For self-averaging quantities an average may not be required. An example is the 
quasi-one dimensional localisation length of the electrons on a very long bar where 
simulation of a single realisation is sufficient. 

This physical quantity T to be analysed depends on the system size L and a set 
of parameters {wi} 



These latter parameters characterise the distribution function of the potential en- 
ergies and also other system parameters such as the energy E, applied magnetic 
held B, spin-orbit couplings, etc. The extrapolation to the thermodynamic limit is 
performed by assuming that that F obeys a scaling law 



Here, for the sake of simplicity, we assume that T is dimensionless. The hope is that, 
in the thermodynamic limit, only one of the many scaling variables (x, </>i, (j>2, ■ • •) 
turns out to be relevant, say \i an d the others {4>i} irrelevant. Here, the words 
relevant and irrelevant are used in the technical sense that the exponent of the 
relevant scaling variable is positive v > and the exponents of the irrelevant scaling 
variables are negative yi < 0. This ad hoc assumption has, of course, to be verified 
during the numerical analysis. 



r = T({w l }),L). 



(1.10) 



(i.ii) 
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For very large systems the contribution of the irrelevant scaling variables can be 
neglected and we obtain a one parameter scaling law 

r = /(£/£), (i.i2) 

with a correlation length, 

s-ixr. (i.i3) 

that depends on the parameters {w,}. This limit is rarely reached in numerical 
simulations and we are forced to deal with the corrections to this one parameter 
scaling behaviour due to the irrelevant scaling variables. (Below we shall refer 
rather loosely to "corrections to scaling" ; strictly speaking we mean corrections to 
one parameter scaling.) 

In practice, we need to simulate not too small systems such that consideration 
of at most one irrelevant scaling variable is sufficient. In this case, the scaling form 
Eq. (jl.lip reduces to 

r = F(xi 1/v ,0L"). (1-14) 

We then fit numerical data for the region close to the phase transition by Taylor 
expanding the scaling function and the scaling variables, and performing a non- 
linear least squares fit. It is important to control the errors in this fitting procedure 
carefully and to specify the precision of all numerical estimates, if the results are to 
be scientifically meaningful. For details we refer the reader to the article by Slevin 
and Ohtsuki. 20 

Such finite size scaling analyses have been used successfully to analyse the An- 
derson transition in three dimensional systems in various symmetry classes, 20,26-28 
the Anderson transition in two dimensional systems with spin-orbit coupling 21 ' 22 
and the plateau transition in the integer quantum Hall effect. 29 

1.3.2. Quasi-one dimensional localisation length 

The next question is which physical quantity to use in the finite size scaling analysis. 
It must be sensitive to the nature, localised or extended, of the eigenstates. (This 
rules out the average of the density of states, for example.) It should also be easily 
determined numerically with a high precision. There are several possibilities. One 
is the localisation length of electrons on a very long bar. Another possibility is the 
level spacing distribution. 30,31 Yet another possibility is the Landauer conductance 
of a hypercube. 26,32 In this section we discuss the first of these possibilities in detail. 

Consider a very long (^-dimensional bar with linear cross-section L. This is a 
quasi-one dimensional system in which all states, irrespective of the values of the 
parameters Wi are known to be exponentially localised with a quasi-one dimensional 
localisation length X(L;wi,w 2 , ...). Using this quasi-one dimensional localisation 
length we define a dimensionless quantity, sometimes called the MacKinnon-Kramer 
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parameter, 

A(L;wi,w 2 ,...) = • (1-15) 

In practice, the error analysis of the simulation is simplified by working directly 
with the inverse of the MacKinnon-Kramer parameter 

T = A-\ (1.16) 

In the localised phase, T increases with L for large enough L, while in the ex- 
tended phase, it decreases. Exactly at the critical point we have scale invariance 
for sufficiently large L 

lim T(L) = const = T c . (1.17) 
1.3.3. The transfer matrix method 

The transfer matrix method is the most efficient way of calculating the quasi-one 
dimensional localisation length. 9 ' 14 The Schrodinger equation for the Anderson 
Hamiltonian on a d-dimensional bar is rewritten as 

V„ : „+ia„ + i = (E — H„)a„ — V„, n _ia„_i . (1-18) 

Here, a ra is the vector consisting of the L^ 1 amplitudes on the lattice sites of the 
cross sectional plane of the bar at n, V rai „ + i is the M x M (M = L d ^ v ) dimensional 
matrix of inter-layer couplings between sites on the cross sections at n and n + 1, 
and H„ is the matrix of intra-layer couplings between sites on the cross section at n. 
Equation (| 1 . 18[) couples the amplitudes of a state at energy E on the cross section 
n + 1 to those at the cross sections n and n — 1. We rewrite (|1.18[) to define the 
2M x 2M transfer matrix, 



V n,n+l( E1 - H ») . - V n,n+l V n,n-l 

1 , o 



(1.19) 



and the transfer matrix product for the whole bar of length TV 

N 

q n = n t « • ( L2 °) 



71=1 

With this, we write 



faN+A Q N ( ai ). (1.21) 



1 I ^6J\ I 

As a consequence of Oseledec's theorem, 33-36 the eigenvalues Ai of the matrix, 

n = ln(Q A rQ] v ) , (1.22) 

obey the following limit 

7i = lim — . (1.23) 
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Here, i indexes the 2M eigenvalues of J7. The values on the left hand side are called 
Lyapunov exponents. They occur in pairs of opposite sign. The smallest positive 
Lyapunov exponent is the inverse of the quasi-one dimensional localisation length, 
i.e. 

1m = \, (1-24) 

where we have assumed that the exponents are labelled in decreasing order. 

Some typical high precision numerical data for the Anderson model in three 
dimensions obtained using the transfer matrix method are shown in Figure ITTTl For 
weak disorder Y decreases, which indicates that in the three dimensional limit the 
system is in the metallic phase. For strong disorder Y increases, which indicates 
that in the three dimensional limit the system is now in the localised phase. At 
the critical disorder, we see that Y is independent of system size. Note that a 
transient behaviour for small system sizes is clearly resolved, which must be taken 
into account by including corrections to scaling when fitting the numerical data. 




4 8 12 16 20 



L 



Fig. 1.1. Numerical data for the three dimensional Anderson model with box distributed random 
potential, width W = 15 — 18 in steps of 0.1. The precision of the data is 0.1%. The lines are a 
finite size scaling fit that includes corrections to scaling. 
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1.3.4. The correlation length 

In addition to the critical exponent v and the scaling functions, one of the principal 
results of the finite size scaling analysis is the correlation length £. We find in the 
localised regime that 

lim \(L) =£. (1.25) 

L— >oo 

Thus, provided the system is in the localised phase, we can identify £ with the 
localisation length in the infinite <i-dimcnsional system. Note that it is important 
to distinguish the quasi-one dimensional localisation A on a long bar, which is always 
finite, from the localisation length £ in the the infinite d-dimensional system, which 
diverges at the Anderson transition. Equation (| 1 . 25|) applies only in the localised 
phase. 

Physically the localisation length £ describes the exponential decay of the trans- 
mission probability t(E; x, x') of a quantum particle between two sites x and x' in 
an infinite d-dimensional system that is in the localised phase 

\ = - lim M^O>. (1 . 26) 
£ |x-x'|->oo I X — x' I 

Thus, the transmission probability, and hence the diffusion constant, vanish in the 
thermodynamic limit and the system is an insulator. 1 

In the metallic phase, the correlation length is again finite and can be related 
to the resistivity. 

1.4. The critical exponents 
1.4.1. Numerical results 

Most strikingly, although corrections to scaling had not been considered extensively 
at that time, already the first works dealing with the orthogonal symmetry class 
showed that the finite size scaling method was able to confirm the most impor- 
tant result of the scaling theory of localisation: whereas in three dimensions clear 
evidence for the existence of a critical point was found, none was found in two di- 
mensions. 9 ' 10 During subsequent years, the universality of the critical behaviour for 
the orthogonal class was explicitly demonstrated by analysing orthogonal models 
with different disorder distributions. 20 It was also demonstrated that the high pre- 
cision estimates of the critical exponents could also be obtained by analysing the 
finite size scaling of various statistics of the conductance distribution. 26 ' 32 

In addition, the critical behaviours of the other universality classes have been 
extensively studied. As can be seen by reference to Table 1, in a given dimension, 
the values of the exponents in the different symmetry classes differ only by several 
percent. Success in clearly distinguishing the critical exponents for the different 
universality classes is a triumph of the finite size scaling method. This is in sharp 
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v = 1.57 ±0.02 


3D orthogonal symmetry 20 


v = 1.43 ±0.04 


3D unitary symmetry 27 


v = 1.375 ±0.016 


3D symplectic symmetry 28 


v = 2.73 ±0.02 


2D symplectic symmetry 21 


v = 2.593 ±0.006 


integer quantum Hall effect 2 



contrast to other methods of estimating the exponents, in particular, the e expan- 
sion, which have singularly failed to yield precise estimates of the exponents and 
even in some cases predicted values that violate the well established inequality 37,38 

2 



v > 



(1.27) 



1.4.2. Remarks concerning experiments 

Measurement of the conductivity at finite temperature on the metallic side of the 
transition and extrapolation to zero temperature permits an estimate of s. Mea- 
surement of the temperature dependence of the conductivity on the insulating side 
of the transition and fitting to the theory of variable range hopping 39,40 permits an 
estimate of v. 

An alternative approach, called finite temperature scaling, 41 is to fit finite tem- 
perature conductivity data on both sides of the transition to 

o-(T) = T s l zv f( X /T 1/zv ). (1.28) 

This permits estimates of s and the product zv. Here, x is the relevant scaling vari- 
able, which is a function of the parameter used to drive the transition. For example, 
for a transition driven by varying the carrier concentration, we can approximate 

X«^), d-29) 
n c 

for doping concentrations n sufficiently close to the critical concentration n c . The 
exponent z, which is called the dynamical exponent, describes the divergence of the 
phase coherence length as the temperature tends to zero 



(1.30) 



Fitting the temperature dependence of the conductivity precisely at the critical 
point, and assuming the validity of Wegner's scaling law Eq. (|1.4j) . permits an esti- 
mate of z. In quantum Hall effect experiments, z has been estimated by exploiting 
the fact that a crossover in the temperature dependence can be observed in very 
small systems when the phase coherence length becomes comparable to the systems 
size. 

The most recent experiments on doped semiconductors 13 have yielded values of 
s and v in the range between 1 and 1.2 that are consistent with Wegner's scaling 
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law Eg. dl.4p . However, there is a clear deviation of the values of v from those in 
Table 1. The most recent experimental estimate of the critical exponent for the 
plateau transition in the integer quantum Hall effect is v = 2.38 ± 0.06. 42 Again 
this differs from the numerical estimate given in Table 1. 

The limitations of models of non-interacting electrons as a description of the 
critical behaviour of the Anderson transition in electronic systems is clearly seen 
in the disagreement between the predicted and measured values of the dynamical 
exponent z. Whereas models of non- interacting electrons predict z — d 43 where d is 
the dimensionality, the experimentally observed value is often smaller. Itoh et al. 13 
found z sa 3 in vanishing magnetic field, which agrees with non-interacting theory, 
but z ~ 2 in applied magnetic field, which does not. For the plateau transition Li 
et al. 42 found z w 1 , which again disagrees with non-interacting theory. 

The advent of experiments with cold atomic gases, 44 and also with ultrasound 
in random elastic media, 45 have allowed Anderson localisation and the Anderson 
transition to be measured in systems that can be reasonably described as non- 
interacting. In particular, Chabe et al. 44 recently measured the critical behaviour 
of the Anderson transition in a quasi-periodic kicked rotor that was realised in a cold 
gas of cesium atoms. For this system, which is in the three dimensional orthogonal 
universality class, Chabe et al. found v = 1.4 ± 0.3; a result that is consistent with 
the numerical estimate in Table 1. 

1.5. Conclusions 

The finite size scaling method combined with high precision numerical simulations 
has permitted the successful verification of the fundamental assumptions underlying 
the scaling theory of localisation and provided high precision estimates of the critical 
exponents. The advent of cold atomic gasses has permitted the experimental ob- 
servation of the Anderson transition in a system that can be reasonably described 
as non-interacting. Describing the critical behaviour observed at the Anderson 
transition in electronic systems remains a challenge and would seem to require the 
development of numerically tractable models that include the long range Coulomb 
interaction between the electrons. 
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